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Abstract 



Based on the non-relativistic QCD factorization formalism, we calculate the bottomonium 
ground state, 775, inclusive charm decays at the leading order in the strong coupling constant 
as and quarkonium internal relative velocity v. The inclusive charm pair production in rji^ 
decay is mainly realized through r]i, ^ c c g process, where the charm and anti-charm quarks 
then dominantly hadronize into charmed hadrons. The momentum distributions of the final 
states are presented. In this work, we also calculate the J/il^ inclusive production rate in the 
r]b decays, where the color-octet contribution is found to be very important. We expect this 
study may shed some light on finding r/^ or knowing more about its nature. 



PACS number(s): 12.38.Bx, 12.39.Hg, 13.20.Gd 
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I. INTRODUCTION 



In high energy research, heavy quarkonium physics is one of the most interesting 
fields and it plays an important role in understanding the hadron configuration and the 
microcosmos. Theoretically, due to the non-relativistic nature of heavy quarknoia, it 
is convenient to research their properties in the framework of non-relativistic potential 
model and non-relativistic QCD []]], or other theories which work well in the non- 
relativistic limit. Experimentally, heavy quarkonia have relatively high production 
rate in both electronic and hadronic collisions and the vector members can be easily 
seen through their bi-lepton decays. Recently, some new resonances have been observed 
in the charmonium energy region 2j, which enrich the heavy quarkonium spectroscopy 
and make heavy- quarkonium physics more interesting. 

After the spin-triplet bottomonium T was discovered three decades ago, its pseu- 
doscalar partner 77^ had been looked for in various experiments. Unfortunately, there 
is still no conclusive evidence that this elusive particle has been found. As a solid 
prediction from the quark model, the existence of rih is indubitable, but its mass and 
decay channels remain undetermined experimentally at the present time. To search 
?7b, several experiments have been conducted both in e~^e~ collisions at the CLEO and 
the LEP, and hadronic collisions at the Fermilab Tevtron. The advantage of e~^e~ 
collisions is its clear background, which is impaired by the fact that production rates 
for spin-singlet states are generally small. Based on the 2.4 fb~^ data taken at the 
T(2S) and T{3S) resonances, CLEO has searched distinctive single photons from hin- 
dered Ml transitions of T{2S) and T{3S) to 77^7, and also from the cascade decay 
T{3S) hfjTT^, /ifcTT+TT^ followed by El transition hf, r]f,'~f, but no signals have been 
found j^. In the experiments at LEP II, the riiy{lS) is expected to be produced in 
two-photon process. The ALPHA Collaborations analyzed the 77 interaction data, 
but found no evident signal in the four- and six-charged-particle final states [4]. The 
results of the L3 Collaboration and the DELPHI Collaboration were also negative and 
the uppe. of the p.oduct= r„ X Br(„) we.e .et fl Q. In comparison to e" 

experiments, the hadronic collision at the Fermilab Tevatron gives large rjf, production 
rate, but due to the complicated hadronic interaction background, the search for rjfy 
there is also hard. Using the full 1992-96 (Run I) data, the CDF Collaboration searched 
the rjb through its exclusive decay to double J/ip and found some oblique evidences. 
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but far from conclusive [8!]. Right now, further efforts is being pursued in the Run II 
data there. 

Considering the situation of rjiy experiments, theoretical research on its properties 
is still necessary, such as its mass, the production cross-sections at different colliders, 
and its various decay channels. Among all properties, the rjfy mass is believed to be 
predicted by potential model, effective theory, and Lattice calculation without much 
ambiguity, which is very important for experimental observation. Recent theoretical 
work fixes the T — rjf, mass splitting in the range of 40 — 60 MeV 



Ref. 



E[lo|,[lll,Q- In 

13|, Braaten, Fleming and Leibovich calculated the rjh production cross section 
at the Fermilab Tevatron in the framework of NRQCD and evaluated the branching 
ratio of the decay rj^ J/ipJ/ip. They suggested that rjh should be observable through 
this decay because of its large branching ratio of 7 x lO"^"*^^. In Ref. Maltoni 
and Polosa also evaluated the observation potential for rjh at the Tevatron, but they 
found and suggested that the decay r]b D*D^*^ might be the most optimistic channel 



to observe rjh signal. In Ref. [15|, the relativistic correction to the rji, — * J/ipJ/ip 
decay process was calculated and a much smaller branching ratio in comparison to 
Ref. [l3'|| was obtained. The author also discussed the r]b D*D^*^ process and got 



a smaller rate than what in Ref. 



141 ]. However, Santorelli finds that the final state 



interaction may enhance r]h J/ipJ/ip decay width by about two orders of magnitude 
[igI . Instead of rjb decays into hadronic final states, Hao et al. calculated the branching 
ratio of the 775 radiative decay process, i.e. rjh — > ''yJ/ip [17]. They claimed that this 
channel is also a hopeful one in the r/;, hunting. 

Comparing to the exclusive process, the inclusive process has a large branching ratio, 
nevertheless normally also has large uncertainties in pinning down the parent particle. 
Fortunately, the final state distributions of experimental observables are helpful to the 
inclusive process for the aim. For instance, recently, the inclusive charm production in 
the Xb and T decays have been studied in [l^ and [3l respectively. 

At the leading order in as and non-relativistic expansion, the ?7f, inclusive open 
charm decay happens through b b g g* followed hj g* —>■ c c, where the initial 
bb is configured in color- singlet. Based on the result of b b ^ g c c, we can roughly 
estimate the branching fraction of b b —>■ X+ charmed hadrons. We calculate this 
process in the framework of NRQCD factorization formalism, in which the un-calculable 
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nonperturbative effects are represented by the matrix elements of NRQCD operators. 
According to NRQCD, the color-octet configurations appear as higher order Fock states 
in rjb decays, which are suppressed by orders of the small magnitude in relative velocity 
f^. Hence, for a leading order calculation, we can safely treat the b b pair inside the 
rib to be in color-singlet. However, for the J/ip production in rib decays, although the 
higher order Fock state, the (ccjl^sf^), is suppressed by v'^ relative to the leading 
color-singlet configuration, it may be compensated by a factor of ag in regarding to 
color-singlet process. 

This paper is organized as follows: the inclusive charm and charmed hadron produc- 
tion in rib decay is evaluated at leading order in Section [ITl In Section UTTI calculation of 
the process rib — > J/tp + X will be presented. The last section is remained for summary. 



II. THE CHARM QUARK PRODUCTION IN % DECAYS 



In this section we calculate the charm quark production in the inclusive rib decays. 
At the leading order in v, according to the NRQCD factorization formulism, its open 
charm decay width takes the form 

T[vb^c-c + X]=C'i (2.1) 

where {0^''{^So)) is a NRQCD matrix element, which represents the long-distance 
effect and gives the probability for finding the heavy quark and antiquark in specific 
configuration within the meson, and can be evaluated by nonperturbative method such 
as lattice simulation. The dimensionless short-distance coefficient can be calculated 
in pQCD. The dominant source of comes from the decay of a color-singlet bb{^So) 
pair into g g*, followed hj g* cc . 

We can calculate the process in the nonrelativistic limit, in which the 66(^S'o) pair 
can be taken as no relative momentum within rib, i.e., Pb = Pb = Q/2, where Q is 
the momentum of the rib. In this situation, for the bb pair to form rib, when it is in a 
color-singlet state, one can replace the product of the Dirac spinors for b and b in the 
initial state with the projector: 

uiPB)v{p-b) ^(^ + M,Jz75 X l-^R^^{0)j ® (^-^^ , (2.2) 



where Nc = 3, and Ic stands for the unit color matrix. M^^^ is the masses of 7]^. At 
leading order in non-relativistic expansion, it could be understood that M^^^ 2 rub- 
The nonperturbative parameters, R{0)rf^ are color-single radial wave functions at the 
origin for rjb, which can be either reached from phenomenological potential models or 
directly extracted from experiments. The relation between the -R(0)^j and {Oi''{^Si)) 
reads (OTi'Si)) = {Nj2n)\R{0\f{l + O{v')). 

One can then get the partial decay width straightforwardly for the process bb{^So) 
c(pi) + c(p2) + gik). That is, 

drfbCSo) ^ ccg] = ^ Yl 1^^*-!' , (2.3) 

ccg 

where M^j^ is the amplitude for the process and the $3 represents the three-body phase 
space, which shows 

In the initial state rest frame, after integrating over the variables which are independent 
of the amplitude , the phase space integration can be further simplified as 

Here, 5*13 = pi ■ k . 

In the numerical calculation, we take mf, = 4.65 ± 0.15 GeV, = 1.50 ± 0.05 GeV, 
and as{mb) = 0.22. The magnitude of the radial wave function at the origin -R(O) of rjb 
equals approximately to that of its spin triplet partner T, which can be determined from 
;he experimental data on the decay width of r(T e+e") = (1.340±0.018) x lO^^GeV 
3- That is: \R{0\f = \R{0)r\'^ = 4.89 ± 0.07 GeVl With these inputs and by 
varying the strong coupling scale from mh/2 to 2mb, for the aim of error estimation, 
we have 

T[bbCSo) ^ ccg] = 190.8™° KeV . (2.6) 

Since at leading order the total decay width of rjb is 

8 

rUvb) ^ rivb - 99) = ^jfrlRv^m' , (2.7) 



6 



3.0x10"* 




FIG. 1: The decay rate variation over momentum fraction yi in the inclusive process % 
c + X, by taking the central values of inputs. 



the branching ratio hence is readily obtained to be 

Br[hh{^SQ) ccg] = 2.61°;^ x 10"^ . (2.8) 

For inclusive decay process, giving out the experiment observable differential dis- 
tribution will be useful. For this purpose, We define two fractions xi = Ei/Ef, and 
Tc = m^/ml, where Ei and E^ stand for the energy of charm quark and bottom quark, 
and in nonrelativitic approximation the Ei, = mi, = The region of variable Xi 

is < Xi < 1. In some cases, instead of Xi it is convenient to use another variable 

which is the momentum of the charm quark divided by its kinematically allowed 



maximum value in rji, decays [18|, ll9| . The relation between these two variables reads 

xi = ^{l-r,)yl + r, , (2.9) 

n - . (2.10) 

The range of j/i is < yi < 1. Figure [T] exhibits the decay rate distribution over the 
momentum fraction yi. 

Because almost all the charm quarks may eventually hadronize into charmed 
hadrons, like D^, D^, D^, and Ac, etc., we schematically show the meson dif- 
ferential distribution in r]b decays in the fragmentation approximation, similar as done 
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in Refs. jisl. It is well-known that the fragmentation function Dc->h{,z) represents 
the probability of a charm quark fragmenting into the charmed hadron h. Here, the 
2; is a Lorentz boost invariant variable, defined as z = In practice calculation, 

we will simply neglect the difference between fragmenting charm quark mass and the 
charmed hadron mass. The z can be reexpressed as: 



(2.11) 



with 



Zh 



a/(1 - r^)yl + re + yi^l - 

1 + v/r^ 



(2.12) 
(2.13) 



i + v^r^ 

Then the momentum distribution of the charmed hadron can be expressed as |18l . Il9l 



dV dzh dzi ^. . .dyi dT 
dyh dyh hu z\ dzi dyi 



(2.14) 



According to the Kartvelishvili-Likhoded-Petrov (KLP) 2l|] fitting for fragmentation 
function 



D,^,{z)=N^z''^il-z) 



(2.15) 



iere, by using the optimal value of ac = 4 for fitted by the Belle Collaboration 



22|. one gets the normalization coefficient A^^ = 8.04 |18j. 



We show in Figure [2] the meson differential production distribution in rji, decays 
in the fragmentation approximation. For other charmed hadrons, the corresponding 
distributions can be obtained similarly. 

It should be mentioned that since the NRQCD factorization breaks down in the 
yi —>■ 1 limit, the velocity and coupling expansions are no more valid in the vicinity 
of the end point. A proper treatment for this endpoint illness is to resum the large 



23 



24 



25 



26 



27 



281] and invoke the shape function 29|. For 



logarithms of log(l — yi] 
these kinds of content readers should refer to the related references in the literature; 
and in this situation, the final state distribution results at the endpoint in this work 
should not be taken seriously. Roughly speaking, when yi < 0.7 the endpoint effects 
become weak and our predictions turn to be robust js^. Fortunately, for total decay 



8 



8.0x10"'' 




FIG. 2: The decay rate variation over momentum fraction in the inclusive process r]i, 
+ X with central values of inputs. 



widths, the endpoint influence is minor and the predictions are quite reUable. For the 
charmed hadron production, apart from the upper point problem, the lower limit also 
poses a problem for the fragmentation approximation, although fl2.15p only comes from 



the phenomenological fitting, 
taken seriously as well. In Ref. 



'herefore, the lower endpoint prediction should not be 



18| , there are detailed discussions about the validity of 



the fragmentation approximation. 



III. % INCLUSIVE DECAY TO J/i/- + X 



In this section, we present the calculation of J/ip inclusive production and its mo- 
mentum distribution in the 77;, decays, as shown in Figure [31 As mentioned in above, at 
the leading order in v only the color-singlet hh(^So) contribution in the initial state is 
the necessary ingredient to be considered, since the higher Fock state contribution can 
not get big enhancement from kinematic or dynamic reasons. However, for the final 
state J/ip both the color-singlet and color-octet effects should be included, because the 
latter one can get one as and kinematic compensation in comparison with the color- 
singlet process. Since it is the hh pair annihilation in the initial state and the creation 
of the cc pair, as well as the emission of gluon or quarks, takes place at the hard scales 
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Vb 



(2) 



FIG. 3: Lowest-order diagrams that contribute to the inclusive process: r]b — > J/ipX. The 
J/^ in Diagrams (1) is a color-single state , while the one in (2) is a color-octet state. 



set by the heavy quark masses, it is legitimate to tackle this semi- inclusive process in 
the pQCD framework. 

The NRQCD formalism allows the systematic calculation of inclusive cross sections 
for quarkonium decays in perturbation QCD to any order in and f ^ , where v is the 
typical relative velocity of the heavy quark inside the quarkonium. Using the NRQCD 
velocity scaling rules l|, we know the color-singlet process, the Figure is at order 
of a^v^, since the matrix element {O^i^Si)) is of order v^. Whereas, the color-octet 
process, the Figure [3](2), is at order a^v'^. Here, v denotes the heavy quark relative 
velocity in Charmonium system. Potential model calculation indicates that the average 
value of v"^ is about 0.3. And, the QCD coupling constant as{mi)) ^ 0.22. In addition 
of the as compensation, the color-octet mechanism may also enhanced by the single 
gluon propagator. Our following calculation really shows that the color-octet process 
should be included in this calculation. 

The decay width can be formulated as 



(3.1) 



where, the Ai and A2 are perturbative calculable short- distance coefficients. We first 
calculate the color-single coefficient. It is customary to start with the parton process; 
here it is b{pb) b{pi) c{pc) c{pc) +c{ki) +0(^2)- Then project the matrix element onto 
corresponding color-singlet configurations. For the bb in initial state, it is the same as 
in the last section. For the J/ip production, the color-singlet projector is 



V{pc) U{pc 



1 



2V47r 



-Rj/i/'(o) 



, (3.2) 
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where e^^^ is the J/ip polarization vector satisfying ej/^(A) ■ ej^^(A') = —5^^' and 
P ■ ^j/ip = 0. -Rj/^(0) is the radial wave function at the origin, which is also valued 
through J/ip leptonic decay width. Combining all these together, one can easily get 
the rjb to J/ip + c + c decay amplitude for the color-single case, i.e. 

X Tr[{Q + M^j757^((fc2 - ^i) ■ 7 + 

1 11 

X -T-, ^-n^ TTTT X 



(^2 - h? - Ml {h + P/2)2 {h + P/2)2 
X u{h)r + Mj/^)Yv{k2) , (3.3) 

where Ci is the corresponding color factor . ki and k2 are the momenta carried by the 
external charm quark and anti-quark, respectively. 

Next, we present the calculation for color-octet process. At the parton level it is 
bipb) KPb) ~^ c{pc) c{pc) + g{k) process followed by projecting the c c spinors onto the 
color-octet configuration, '^Sf\ while keeping on configuring the initial 6 6 in color- 
singlet. The color-octet projector is 

v{p,)u{p,) — ^ fj,^{r + Mj,^) X (^^=R^j^^{Q)\(^V2^ . (3.4) 



Here, T^j denotes the SU(3) generator. And, we introduce another phenomenologi- 
cal parameter i?jy^(0), which stands for the color-octet nonperturbative effect. The 
relation between Rj^^{0) and the NRQCD matrix element {Og^^ {^Si)) is defined as: 

{OrCSr)) = ^-^m/40)r (3.5) 



From the fitted value of {Os^^{^Si)) ^ 1.5 x IQ-^GeV^ |3l|, we have \R 



j/i, 



0.102 GeV'^/^. Then the decay amplitude for color-octet case is 

«^ ^ ,^r,,(0)i?^,/,,,(0) 



X Tr[{Q + M,j75 /;(A2)((Q/2 - fc) ■ 7 + M,,/2)7l 

xrr[/j/^(Ai)(/' + Mj/^)7j , (3.6) 



where Cg is the color factor, k is the momentum carried by the external gluon, and 
is the gluon polarization satisfying ■ e^, = 0. 
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With the matrix elements Mj^^ and Mfj^ we can immediately get the rji, J/ip + X 
decay width. The analytical result for it is a bit lengthy, and will not be presented 
here. For our numerical estimation, the non-relativistic limit is also enforced for the 
charmonium. That is we take the Mj/^ ^ 2mc approximation. From T{J /ip — > e^e~) = 
(5.55±0.14) X lO-^GeV [201, we get |i?j/^(0)|2 = 0.527±0.013 GeV^. Then, the decay 



width for the concerned process reads, 

r(r/, ^ J /^color-singlet + X) = Q.U+J,^^ KeV , (3.7) 
r(r76 ^ J /^color-octet + X) = 2.16 KcV , (3.8) 

^totai{^h ^ J/i^ + X) = 2.29 KeV . (3.9) 

That means that the r/;, J/ip + X process has a branching ratio of 3.23 x 10~^ 
or so in the r^b decays. Here, the uncertainty estimate of the color-singlet process is 
performed in the same way as in the preceding section. Whereas, considering of the 
large uncertainties remaining in the color-octet matrix element fitting, we carry the 
numerical calculation for color-octet process by only taking the central values of the 
inputs. 

Like in section II, to give out the differential decay width we define three new 
variables, the X2, rjj^, and 2/2 , as 



X2 = Ejl^/Ey, , (3.10) 



Tji^ = M],^/ml , (3.11) 



y^ = xr-t^- (3.12) 
Then we can express the partial decay width of rji, — > J / ip color -singlet + X as 
dT[rn, J/^coior —single 

+ ^] = ^ E \^ltr?d^, (3.13) 



In analogy to what is performed in the last section, we get the momentum distri- 
bution dVirji, — > J /ip color -single + X) / dy2, as showen in Figure HI For the the process 

Vb ~^ J /i^ color -octet + X, 

dr[7], ^ J/Aolor-octet + X] = ^ ^ \M^tri>^,, X2)\'d<^2 • (3.14) 

Ai,A2 
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FIG. 4: The J /ijj momentum distribution in the inclusive process % J /i^ color -singlet + X. 



Since this is a two-body decay process, the J/ip momentum distribution dT{rii, — >■ 

J li^ color-octet + X)/dy2 is ouly a delta function peaked at ?/2 = W 4^2 (^2 -4m^^^) = 
Again, for the to J/V^ inclusive decay distribution, one should pay attention to the 



endpoint problem 32] • In particular for the color-octet contribution, the r^f, two-body 



decay at leading order resulting in a delta function distribution, which is smeared out 



by the non-perturbative effects and resulting in a shape function [33|. 



The numerical result shows that the branching ratio for color-octet process is larger 
than the one for color-singlet process by about an order, which offers an opportunity 
to check the existence of color-octet mechanism experimentally. Considering of the 
uncertainties existing in the magnitude of color-octet matrix element, the numerical 
difference between these two processes might shrink down, nevertheless, they give a 
distinctively different momentum distribution, which may also help experimenters to 
distinguish them in the future experiment. 



IV. SUMMARY 



We have studied in the framework of NRQCD the inclusive charm production in 
the decay of the pseudoscalar bottomnium state rjb. We find that it gives a quite 
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large branching fraction. Since the produced charm quarks will dominantly evolve 
into charmed hadrons, by employing the KLP fragmentation function, we give out the 
momentum distribution of D'^, as an example. 

We have also calculated the decay width and the momentum distribution of the 
inclusive J/ip production in the rjb decay. We find that in this case the color-octet 
process should be taken into consideration. However, this two different J/ip production 
schemes have obviously different momentum distributions. This is a distinct character 
of this process, which will be helpful for future experiment to investigate the rjb and to 
the study J/ijj production. 

In all, to investigate the elusive 77^ is still an interesting task for both theory and 
experiment. Our explicit calculation shows that rjb inclusive decays to charm pair (in 
experiment the charmed hadron pair) and J/ip have quite large branching fractions. 
These processes can be helpful for people to hunt for the rjb at the Fermilab Tevatron, 

or LHC, where copious rjh are expected. 
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